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We study quantum optics with the atoms coupled to the quantized electromagnetic (EM) field in a non- 
inertial reference frame by making use of quantum field theory in curved spacetime. We rigorously establish 
the microscopic model for a two-level atom interacting with the quantized EM field in a rotating ring cavity by 
deriving a Jaynes-Cummings (JC) type Hamiltonian. Due to the two fold degeneracy of the ring cavity modes, 
the Rabi splitting exhibits three rather than two resonant frequency peaks. We find that the heights of the two 
side peaks show a sensitive linear dependence on the rotating velocity. This high sensitivity can be utilized to 
detect the angular velocity of the whole system. 
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Introduction - The interference of of two light beams in a 
rotating ring can be utilized to measure the rotating velocity. 
This is well known as the Sagnac effect, which bases some 
optical gyroscope schemes [1-6]. For the Sagnac effect in the 
medium with linear dispersion, there have been a lot of so¬ 
phisticated studies based on classical optics. If we want to 
make the optical gyroscopes to an extremely high precision, it 
is necessary to consider the quantum fluctuations in this rotat¬ 
ing optical system. We notice that a rigorous quantum theory 
about the microscopic model about the interaction between 
the atoms and the quantized EM held in a rotating reference 
frame is still not well established, but it is obviously essential 
for the study of quantum and nonlinear effects in a rotating 
optical system. 

In this letter, we ascribe the effects of rotation to the 
“curved” spacetime metric according to the generic principles 
in relativity. Starting from the classical Lagrangians of the 
EM held and a charged particle based on the the principle of 
least action, we obtain the covariant motion equations in the 
rotating non-inertial reference frame (NIRE). The variation is 
carried out in the rotating reference frame, and all the non- 
inertial physical effects rooted in the rotation are included in 
the “curved” spacetime metric. Then we obtain the quantized 
Hamiltonian for quantum optics in NIRE through the canoni¬ 
cal quantization. 

We derive a microscopic model of an atom interacting with 
the quantized EM held in a rotating reference frame, which 
then gives the Jaynes-Cummings (JC) model for a rotating 
ring cavity coupled with a two-level atom. Eor a ring cavity 
in the inertial rest frame, the clockwise (CW) and counter¬ 
clockwise (CCW) propagating optical modes are always ex¬ 
actly degenerated. Thus, in the JC model of a ring cavity, 
the atom couples with the two degenerate modes simultane¬ 
ously. Due to the existence of two optical modes, the Rabi 
splitting of this system exhibits three rather than two resonant 
frequency peaks. More importantly, we And that the rotation 
of the ring would induce a detuning between the original de¬ 
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Figure 1: (Color online) Schematic setup. A two-level atom is fixed 
in the rotating ring, and external fibers are used for driving and prob¬ 
ing. Two photon detectors (PD) are used to measure the photon cur¬ 
rent. 


generated modes in the non-inertial frame, and the heights of 
the two side peaks would change with the rotating speed due 
to this rotation induced detuning. At lower speed, the heights 
of the side peaks depend linearly on the rotating velocity. This 
sensitivity can be utilized to detect the rotation of the whole 
system, which can be regarded as a quantum Sagnac effect. 

EM field in the rotating reference frame - We consider 
the EM field rotating around the z-axial in the CCW direction 
with angular speed H. Physics laws must have the same co¬ 
variant mathematical form in all reference frames (including 
the non-inertial rotating frame that we are studying). Thus, 
in the rotating reference frame, the Lagrangian density of the 
EM field is 

C = ( 1 ) 

4/ro 

Here po is the magnetic constant, := = 

df^Ai, — dyA^ and F^'' := F^p. In the above def¬ 

initions, dfi := djdx^, where := {ct,x,y,z) is the 4- 
dimensional coordinate, and A^ := {—(fjc, A/ (i = 1, 2, 3) 
is the electromagnetic 4-potential. is the covariant deriva- 
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tion defined by 


=df.A,-Tl,Ar. (2) 


Here is the Christoffel symbol and it can be calculated 
from the spacetime metric [7, 8]. All the physical effects 
due to the rotation are included in the “curved” spacetime met¬ 
ric (see Ref. [2, 9] or Appendix A), i.e., 


= 


-1 

C 

Qx 


1 - ^ 
Vrxy 


Q.X 

xy 

-pr- 


1 - 


a^x^ 

-pp- 


(3) 


From the variation of the action 
6S = . 


(5 y £ • ~ ^ J 


we obtain the Euler-Lagrangian equation as 


dC 


d{dy,A^ 


1 


dC 


dA^ 


(4) 


= 0 , 


where g := det g^". Here C := is a functional of 

and dfiAy. This is the covariant form of the Maxwell equa¬ 
tion, which applies to all general reference frames (both iner¬ 
tial and non-inertial ones) [7]. Substituting the spacetime met¬ 
ric Eq. (3) into the above covariant equation, and taking the 
Coulomb gauge V-A = 0, Aq = 0, we obtain the d’Alembert 
equation = 0 as 

[-dl + • Vao - {^R ■ Vf]A, = 0. (5) 

Here := c~^dldt, vr := vr/c, and vr := fi x r is the 
linear velocity. 

Next we reduce the problem into a quasi-ID ring configu¬ 
ration. We assume that A(r) is homogenous in the transverse 
direction, and we have V ~ es9s [9], where bg is the direction 
along the ring. Then we have 

[-So-b 2{;fl5g5o-b (1 - u|)9g]Aj = 0. (6) 


Quantization of the EM field - Under the Coulomb gauge 
we used above, the canonical momentums are = 0 and 

dC 

where Eq is the electric constant. The Hamiltonian density 
Tf = Ai — £ is obtained as 

7f = i(EoE2-b—B2)-lEo(vflxB)2. (10) 

Z flQ Z 

Restricted in the quasi-ID ring configuration, respectively 
they become 


£* = Eo(StA - VRdsA)i, 

7f = i£o(atA)2 + i^(SgA)2. (11) 

Z Zgo 

Notice that since V ~ e^Ss, the Coulomb gauge V • A = 0 
leads to A • §8 =0. Namely, A(s, f) only has two transverse 
directions, and so does the canonical momentum £{s, t). 

We apply the following canonical quantization condition 

[Axis, t), S'^is', f)] = ihSlSis - s') ■ (12) 

[AA(s,f), A^(s',f)] = [£^(s,f), £‘^(s',<)] = 0. 

Here X,a = 1,2 mean the two transverse directions, and S is 
the cross-sectional area of the ring. Since we have reduced the 
problem into 1-dimension the above quantization condition is 
consistent with the Coulomb gauge condition V • A = 0 au¬ 
tomatically. Now we write down the field operator A(s,f) 
as 


A(s,<) =^efcAZfc[afcAe*'=*-“'=* + 4;,e-*“'=‘]. (13) 

k,X 

With the help of the above canonical quantization condition 
(12), we can prove the following bosonic commutation rela¬ 
tions (see Appendix B), 


Here vr = vr/c and vr := |v/{| = flR is the linear speed. 
The general solution of the above equation is 

A(s,t) = (7) 

fc,A 

where k = 27rn/L, n G Z, and L is the length of the ring [14]; 
ekx is the polarization directions in the transverse section, and 
Zfc is a normalization constant. The eigenvalue equation of 
Eq. (6), ujI -b 2vRkuJk — (c^ — vj^)k^ = 0, gives rise to the 
following dispersion relation 

Wfc = ±ck - VRk. (8) 


[OfcAi dqfj] — Ago-] — 

[ofcA, d\^] = SkqSx. ■ ■ (14) 

Thus, the normalization constant is taken as = 

^/h/2eoVc\k\, where U = L • S is the effective volume of 
the ring, so that [a^A, ajcr] = dkqSxcr- Then we obtain the 
quantized Hamiltonian of the EM field as 

Hem= [ dVU^^nwkid\^dkx + ^). (15) 

fc,A 


Therefore, the rotation makes the dispersion relation 
anisotropic, i.e., for the two modes = jfcj and fc_ = —j/cj, 
their frequencies ujk^ no longer equal. 


This Hamiltonian has the same form as that of the EM field 
in the inertial frame, but the dispersion relation Wfc is changed 
due to the rotation [see Eq. (8)]. 
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JC-model in the rotating frame - Next we study the mo¬ 
tion of a charged particle in the EM field in the rotating refer¬ 
ence frame. To this end, we start with the invariant Lagrangian 
of a charged particle in the effective curved spacetime [7, 8 ] 


. dx^ dx ’', 1 . dx^ 

L=-n,c[-s,,——\n+eA,—, 


(16) 


where r is the proper time. The Hamiltonian description is ob¬ 
tained from the following variation (see Ref. [ 8 ] or Appendix 
C) 


5S = 




S / Ldt 


Here := dx^/dt, : 
and 


= S j [P,v^ - He]dt. (17) 
dXfj,/dt = L := L ■ T'^ 


r 


dt c 

dr \/—vov^ — UiU® 


(18) 


As a conservative system, the Hamiltonian He is the func¬ 
tional of x^(t) and Pi(t) (i = 1,2,3), and 


dL 

Pi = TTT = Tmvi + eAi := pi + eAi, 
ov^ 

He = Piv'' — L = —Tmvov^ — cAqv^. 


(19) 


Here pi := TmVi is the mechanical momentum. Replacing Vi 
and u* by the canonical momentum P^ = pi + eAi in Hf., the 
Hamiltonian of the charged particle is obtained as (see Ap¬ 
pendix C) 

He = - g°°{ 9 "^P^P 3 + m'^c'^)]-ev°Ao. 

Substituting in the metric [Eq. (3)], the above Hamilto¬ 
nian becomes 

He = v/(P — eKY(? -I- -|- v/j • (P — eA) -|- ep 

~ ^(P - eA)2 -b • (P - eA) -f ep, (20) 
zm 

where P := (Pi, P 2 , P 3 ). In the quasi-lD ring, A is always 
in the transverse direction and perpendicular to the linear ve¬ 
locity v/{, thus V/j • A = 0. 

We consider an atom with the nucleus fixed at a certain po¬ 
sition inside the ring cavity. Around the nucleus, the electron 
is trapped by the central force potential p ~ r~^. Thus, in the 
expansion of the above equation, Patom := P^/ 2 m + ep + 
VR ■ P is exactly the Hamiltonian of a hydrogen-like atom plus 
a perturbation term vr ■ P, which comes from the non-inertial 
effect of rotation, and the term P • A describes the coupling 
between the atom and the EM field. 

Treating P and ic as operators, we obtain the energy lev¬ 
els of the hydrogen-like atom. Then we focus on two en¬ 
ergy levels |e) and |g) with a dipole transition. Neglecting 
the A^ term, with the help of Eq. (13), the total Hamiltonian 
H = He + Hem reads, 

H = gk\dkx+V*kxdl^]+'^f>^kdl^dkx, 

fc,A fc,A 


where := |e)(e| - |g)(g|, := |e)(g| -b |g)(e| and &y = 

i[a,^ d^]. The coefficients ^ and r]kx are explicitly given as 


( = i(e|vi^ • P|g) = 


Vkx = —{e\ZkP^^(F-ekx)\g) 
m 


( 21 ) 


where p := e(e|x|g) is the dipole moment [10]. Here the 
dipole approximation is applied, and sq is the position of the 
atom. We choose sq = 0 to cancel the phase factors. Com¬ 
paring with the inertial case vr = 0 , we see that fl and gkx 
is unchanged, and the contribution of rotation appears in the 
correction term ^d^ and the dispersion relation Wfc. 

We set the two polarization directions to be parallel and 
vertical to the projection of p in the transversal section re¬ 
spectively. Then the vertical optical modes are decoupled 
with the atom [see Eq. (21)]. When the ring is not too long, 
the frequencies of different eigen modes of the EM field 
are well separated from each other, so we only consider the 
modes nearly resonant with the atom. But we should notice 
that, different from the Eabry-Perot type, the = |A:| and 
k- = —\k\ modes are always nearly degenerate in ring cavity 
(unless k = 0). When vr = 0, they are exactly degenerated 
. That means, the and fc_ modes must be con¬ 
sidered together. Therefore, by omitting the double creation 
and annihilation terms, we establish the JC-model of two-level 
atom in a rotating ring cavity with the JC Hamiltonian 


- kPh ^ ^ .v'I'a ..'I'., 

H + Iia;+d]|_d+ -b /ia;_a]_a_ 

+ gd'^{d++ d-) + g*d~{d\+d}_), ( 22 ) 


where the coupling strength g is 


g = in{p-ekx)[ 


h 

2eoVc\k\ 


1 

2 


(23) 


Here d± is the annihilation operator for the modes k±, and 
their frequencies are u;± = (c=F VR)\k\ := ujo ± A, where we 
define wq := c\k\ as the rest frequency of the ring cavity and 
A := —VRk as the rotation detuning. Eor simplicity, hereafter 
we absorb the phase of g into the operators to make g = g* > 
0 and set ti= 1. 

Rabi splitting in the rotating frame - We consider the 
case that the dipole moment is parallel to the transverse direc¬ 
tion, thus ^ = 0 and we can choose the polarization direction 
to satisfy p • ekx = |p|- In this case, the excitation number 
N := d^ + al,_a+ -b a’Ld- is conserved, i.e., [TV, H] = 0. 
The ground state is |G) = |g, 0,0) and the eigen energy is 
Eq = —H/2. But generally we cannot give an analytical so¬ 
lution for the whole energy spectrum. When ujq = H, we 
can obtain the eigen energy of the first three excited levels 
[Eig. 2(d)], and they are 


En = 


n 


n 

E± = -±A„ 


(24) 
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Figure 2: (Color online) (a, b) The steady photon number n± of 
the two cavity modes. The two-level atom is resonant with the rest 
frequency 12 = wo, and we set wo = 12 = 1 as the unit. The other 
parameters are g = 1 x 10”"^, 7 = 0.5 x 10”“*, S = 0.05 x 10”"^, 
and (a) A = 0 (b) A = —VRk = 1 x 10“®. (c) The height of the 
side peak of n+(S2) at 12 = ^/2g changes with the rotation detuning 
A. (d) Demonstration of the ground state and the first three excited 
states. 


where Ag := + 2g‘^ (see Appendix D or Ref. [11]). 

Therefore, if we use an external input to drive the ring cav¬ 
ity weakly, we predict to see Rabi splitting with three resonant 
peaks [12]. We consider a probing setup as demonstrated in 
Fig. 1. An external driving laser is input to drive the fc+-mode 
of the ring cavity. The photons in the ring cavity can leak into 
the output fiber, and then be probed by photon detectors. 

We use the following master equation to describe the sys¬ 
tem, 

p= i[p,H+Hd{t)]+ ^ ^(2aapa]j-{p,a]jaa}), (25) 

a—+,- 

where Hd{t) = is the driving term. 

Here we consider the case 12 = wo- When the driving strength 
£ is weak, the cavity modes will not be excited to states with 
large photon numbers, and we obtain (see Appendix D) 

n+ ~ —, n_ ~ 16£ (26) 

where 

M ■.=4[g^ - - A)f +, (27) 

F + 87^[2/ + 12^ - 

+ 16122[2p2_H2_^A2]2, 

and 12 := 12 — Wd- 

When 7 —^ 0 , the denominator F(12) has three minima 
around 12 = 0 and 12 = ±y/A‘^ + 2g‘^, which give rise 
to three peaks in n±( 12 ) corresponding to the first three ex¬ 
cited energy Eq. (24). We can also explicitly see that F(12) 


is symmetric for ±12, i.e., F(Cl) = F{—Cl), but M(12) is 
not. Therefore, n_(12) is symmetric for ±12, but n+(12) is 
not. We plot the steady photon number 7 Td-( 12 ) of the cav¬ 
ity modes in Fig. 2. When there is no rotation, vr = 0, the 
steady photon number of the two modes n± ( 12 ) are both sym¬ 
metric with respect to the driving detuning ±12 [Fig. 2(a)]. It 
is worth noticing that when there is a small rotating velocity, 
the heights of the side peaks of n+( 12 ), which correspond to 
the mode a+ being driven, change sensitively with rotating 
velocity [Fig. 2(b)]. 

Since the positions of the side peaks are 12 ~ 

±y/A 2 TV, we plot the height of (12 = \/2g) with re¬ 
spect to the rotation detuning A = —vnk around A ~ 0 
[Fig. 2(c)]. When A is quite small (in this example, we 
have |A/12| = |ffi/c| < 10“®), n+(A) depend linearly on 
A = —VRk around A = 0. From Eqs. (26, 27), we obtain the 
slope of n+(A, 12 = '/ig) around A = 0 as 


cM+(12, A) 


dA 




MV2g£‘^ 

72(72 -F 852)' 


(28) 


From these results, we see that the sensitivity of this height 
change with respect to the rotating velocity can be controlled 
by the coupling strength g and the decay rate 7. A ring cavity 
with high quality promises a sensitive measurement. 

In experiments, the steady photon number n± can be mea¬ 
sured by the average output photon current. In a probe setup as 
shown in Fig. 1, we have a±_ouT = a±.iN ± v/7®± ’ 
average output photon current is (a^ out®±,OUt) = 'yn± 
[13]. This photon current directly characterizes the steady 
photon numbers of the cavity modes, and can be measured 
by the photon detectors. 

Summary - In conclusion, we have generally studied the 
quantum optics for the interaction between light and atom in a 
non-inertial reference frame by the approach of quantum field 
theory in curved spacetime. For the two-level atom interacting 
with the quantized EM field in a rotating ring cavity, our ap¬ 
proach built a microscopic model described by a two-mode JC 
Hamiltonian. Based on this generalized JC model, our study 
predicts that that the heights of the side peaks in the Rabi split¬ 
ting show a sensitive linear dependence on the rotating veloc¬ 
ity at low speed. Therefore, this model can not only be utilized 
for hybrid optical gyroscope design, but also provide the fu¬ 
ture development of quantum gyroscope with a solid physical 
base to consider the effect of quantum fluctuation. 
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where 77 ^“^ = diag{ — 1 , 1 , 1 , 1 } is the metric in the inertial lab 
frame. And we also have 
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Appendix A: spacetime metric 

Considering the system rotating along the z-axial in the 
counter-clockwise direction, we have [2, 9] 



C 


c 


c c 

1 

1 

1 


(A4) 


Appendix B: Canonical quantization of the rotating EM field 

In the canonical quantization of the rotating EM field, it is 
crucial to find a proper orthogonal relation of the field oper¬ 
ator A(s, f). This relation can be obtained with the help of a 
density-flow relation derived from the equation of motion. 
From the equation of A(s, f) 

[-^2 + 2iRdsdo + (1 - = 0, (Bl) 

we obtain the following density-flow relation 

5o[A*27A] = a,[A*:TA], (B2) 

where 

A*VA :=(A* • doA - A ■ doA*) 

-vr{A* -dsA- A-dsA*), 

A*JA :=(1 - fl){A* -dsA-A- d^A*) 

+ vr{A* -d^A-A-doA*). 

Integrating Eq. (B2) over the whole space, we obtain 


t = T, 

a: = A cos Of-f Esin Iff, (Al) 

y = —X sin ilt -\-Y cos Of, 
z = Z. 


Here, X°‘ = [cT, X, Y, Z) is the coordinates in the inertial 
lab frame, while a;“ = (cf, x, y, z) is in the co-rotating frame. 
We obtain the general coordinate transformation matrix as 


[A“/3] 


dx 


1 

Oy/c cos Of sin Of 
—ffxjc — sinOf cosOf 

1 


■ (A2) 


The spacetime metric in the co-rotating frame, as a co¬ 
variant tensor, can be calculated by 


[9^ = [ 


dx^ 

-1 

Oy 

c 

Qx 


■V 


a/3 


dx’' 


_Qy 
1 _ 

-L _9 


Qx 

xy 




1 - 




(A3) 


do(^J dS- ds[A*VA]^ = 0, (B3) 

where f dS is the integral over the cross section and gives a 
constant area S. Then we know that this integral is a conserved 
constant independent of time (but still not determined yet). 

Then we can find some orthogonal relations. For the eigen 
solution 


AfcA(s,f) := efcAafc(s,f) := CfcAZfce*^® (B4) 


we have the following orthogonal relation 



2f|ZfcpL 

c 


(Wfc +VRk) SkqSxa, 


ds[A;,VAl^] = 0. (B5) 


Notice that no matter k > 0 and fc < 0, we both have ujk + 
VRk = c\k\. 

With the help of the above orthogonal relation, for the field 
operator 


A(s,f) = ^efcAZfc[afeAe*'=*-'^'“‘ + 

fc,A 
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we have 


ds[Al^VA] = -2i\Zk?L\k\-akx, 


f 

[ ds[AkxVA] = 2i\Zk\\\k\ 
Jo 


'A\- 


(B6) 


Notice that indeed the above terms in the integrals can be 
written as 

■ {do - VRds)A - A-{do- VRds)&l)^ 

= — i\k\Zf.A\^ 

CEo 

AkxDA = ■ {do - VRds)A - A - {do- VRds)akx 

ak£^ 


ceo 


+ j|A:|afcAA. 


Then, together with the canonical quantization condi¬ 
tions [Ax{s,t), £‘^{s',t)] = iliS'^S{s — s') ■ and 

[Aa(s, t), Acr(s', <)] = t), £'^{s', t)] = 0, we can cal¬ 

culate [dfcA, ajcr] from Eq. (B6) as follows 

[ 2 |Zfc| 2 L|A:|.afcA, 2 |Z,pL|g| • a+J 

= JJ dsds'[Al)^{s,t)'DA{s,t), Aqcr{s',t)'DMs'^t)] 


where H[xi{t), P'{t)] is the functional of independent vari¬ 
ables x®(f) and Pi{t). The above variation gives 


0=(5 


(]t^ \ 

P^-H[x\P]) 


dt 


dt 


(C2) 


=P^6x' 


B 


f.dx^ dH dPi dH 

+ - a?; + 


A and B are the initial and final points of the trajectories of 
x®(f) which are fixed, and thus the first term is zero. Since 6x' 
and SPi are independent variables, the above variation leads 
to the Hamilton equations 


dx^ 

dt 


dH dP, 
Wi' ~dt 


dH 
dx^' 


(C3) 


Here we have utilized the following relation 


■I 


6 / Piv'^dt = 


P5P,+P,5{—)^ 


v'^SPi + Pi — iSx'') 
dt 


dt 

dt 


= PAx' 


/ 


J-x'' ■ SPi — '^Pi ■ Sx^ 
Idt dt 


(C4) 

dt. 


JJ ' C£o 

2 ^L|fc| - 2 


C£o 


- , , Therefore, for a charged particle in the co-rotating EM field, 

Aq\a,A^{s )]^,have 


C£oS 




AOfeg- 


Therefore, we choose the normalization constant Zfc to be 
h , 1 , fi , 1 


Zfc = 


2eoV{ujk + VRk) 2£oEc|fc| 


(B7) 


dL 

Pi = ttt = Tmui -I- eAi := pi -f eA^, 
dv^ 

A 

H = Py-L = Tm{vy + ^) - eAov° (C5) 
= —Tmvov^ — eAov^. 


where V := L• S is the volume of the quasi-1D ring, so that we 
obtain the bosonic commutation relation [dkx,oJqa] — SxaSkq- 
With the same method, we can also check [dkx, dqa] = 



After we have got the normalization constant Z^, it is 
straightforward to verify that the Hamiltonian of the rotating 
EM field is 

Hem = I dVn 

= I dE[ieo(5,A)2 + i^(a.A) 2 ] 

= \f'^k{dl)^dkx + akxdiO- 

k,X 


Appendix C: Hamiltonian description of a charged particle in 
the co-rotating EM field 


The Hamiltonian equation can be derived from the follow¬ 
ing variation. 


■jMr.ij. 


dx^ 


Ldr = S I Ldt = S {Pi—- H)dt = 0, (Cl) 


Here we denote := dxk-/dt,v^ := dx^/dt. And we denote 
Pi := Tmvi as the mechanical momentum. Here we take 
xk" = {ct, xi,X 2 ,X 3 ) and = Pfiyx'', thus we have = c. 

Eurther, we still need to replace E mvo by the canonical mo¬ 
mentum Pi in the Hamiltonian H. We emphasize that here 
we only have got the relation between Pi, pi and Vi from 
Eq. (C5), but we did not have the definition of Pq, po, so we 
cannot just naively replace Emuo as po- To find the relation 
between E mvo and Pi, pi, we should notice the following two 
relations. 


r^m^vy = ym^Vi{g^°Vo + g'^Vj) (C6) 

= 9'°Pi ■ Tmvo + g'^iPj 


t'2 2, i 2 2 

i m ViV = m c 


ViV" 


— VoV° — ViV'^ 


(C7) 


= —m^c^ — Tmvo ■ Tmv^ 

= — (?°°(Emuo)^ — ■ (Emuo) — 


Therefore, we obtain an equation about Emuo, i.e., 
g°°{Tmvo)‘^ + 2g°'’pi ■ (Etouq) = 0, (C8 ) 
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which leads to the solution (Pmuo should be negative, so the (t^/ 2)], and the equation becomes time independent. Then we 
other positive solution is invalid) obtain equations for observable expectations as follows 


+ rri^c^) 

Tmuo = --^55-. 

Therefore, we have 

He = - g°°{9^^PiP3 + m‘^c^)]-ev°Ao, 

where Pi = Pi — eAi. 


Appendix D: Steady photon number 


For the JC-model for a two-level atom in a rotating ring 
cavity 

'' /iff + + 

77 + ^(7^ +/ia;_|_a^a++/ia;_aLa- (Dl) 

-I- p6-+(d+ + a_) + g*a~{d\_ + d^_), 


we can directly verify that the ground state is |g, 0, 0) and the 
eigen energy is Eq = —il/2. 

The excitation number N := + d^a+ + d^_d- is al¬ 

ways conserved, i.e., [N, H] = 0. Thus, the single excitation 
subspace, which is spanned by |e, 0,0), |g, 1,0), |g,0,1), is 
closed, i.e.. 


77|e,0,0) = ^|e,0,0)+5|g,l,0)+p|g,0,l) 

H\g, 1, 0) = 5|e, 0, 0) + (^^ + A)|g, 1, 0) (D2) 

77 |g,0,l)=5|e,0,0) + (^^-A)|g,0,l) 

When n = Wo, the above eigen equations can be diagonal¬ 
ized exactly and we can obtain the eigen energy states in the 
single excitation subspace. The eigen energies are T/q = 17/2 
and E± = f2/2 ± Ag, where Ag := and the 

eigenstates are 


\Eo) [A|e, 0,0) - 5|g, 1,0) + g\g, 0,1)1, (D3) 

Zo L J 

\E±) ~ [^(A ± Ag)\e, 0,0) + i(A ± A,)2|g, 1,0) 

+ l<?ng,0,l) 


where Zq and Z± are normalization constants. 

When we use an external laser to drive the d+ mode, the 
master equation of the system is 

p = i[p,H+Hd{t)]+ '^{2dapd'l-{p,dida}), (D4) 

where is the driving term. We 

make a unitary transformation by exp[zwd(dl|_a+ + dPd- + 


0 = 

-i{dj+ 

- 


0 = 

—i{uj- 

-/-)a- -ig{d ) 


0 = 

-i{ui+ 

-i^)S+-i£{d^)+ig{d' 

-) 

0 = 

—i{uj- 



0 = 

-/(w+ 

— Ft — i^)Z+ — iS{d^) — 

■9, 
* 2 ' 

0 = 

—i{uj- 

-h-^l)Z_-^^-{{d^)Fl) 

0 = 

1 

. 

) + ig{s+ + S-) 


0 = 

(Z+- 

Zl) + (Z_ - Z*_) 



Here we denote Z± := {d^d±), S± := {d^d±), a± := 
(d±) and w± := w± — Wd, If := H — Wd. In this 
equation we have omitted terms of higher orders, like 
((T^dlj_d±), {d^d\d-), so the above linear equations be¬ 
come complete. This approximation requires that the driving 
strength is weak thus the excitation is very low. 

From this set of linear equations, we obtain the solution of 
a± = (d±), and then we obtain n± ~ \a± p as follows 


where 


dS^[MF'^D + o{S^)\ 
{ED -b f ’ 
WS'^g^ ■ E^D 
~ -b 4£’2G')2 ’ 


(D5) 


M ■.=A{g^ - HA _)2 -^ 7 ^ 0 ^ 

D 1=7^ -b 4A^^ 

F -b 472(2/ -b (A+H - g^f -b - g^f] 

-b 16 H 2 [ 2 / _ ^2 + A2]2^ 

G :=16A2 [( 2 / - ^ 2)-2 2p2f)A] -b 7 '‘( 2 / _ ^ 2 ) 

-b 472 [(2/ _ ^2)^2 ^ 2^,2 _ ^2^^2 2 /HA] . 

When the driving strength E is weak, we omit terms of o(£’2) 
in Eq. (D5) and obtain 


n+ 


4£^M 

F 


71- 


16£’2p'^ 

F 


(D6) 


For the case that the atom is resonant with the resonant fre¬ 
quency, i.e., O = Wo, we have w± = H ± A, and we obtain 

M := 4 [/ _ Cl{Cl - A )]2 -b 721)2^ 

F + 872(2/ -b - 252^2 + A 2 h 2 ] 

-bl6H2(2/-H2_^ A2]2 


This is what we have shown in the text. 

The above master equation can be also solved numerically 
by setting a cutoff on the Hilbert dimension of the cavity 
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modes. We find that when the driving strength £ is weak, well accordance with the numerical calculations, 
the above analytical results of the steady photon number show 


